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’JP ?£? -ATH0 SPHERE  DENSITIES  DURING  AN  INTENSE  MAGNETIC  STORM, 
FROM  THE  LOGACS  EXPERIMENT 


by 

R.  P..  Allan 
G.  E.  Cook 

SUMMARY 

One  of  the  most  intense  magnetic  storms  of  recent  years  occurred  on 
25-26  May  1967  and  coincided  aLt*i  the  flight  of  the  LOGACS  (Lov-G  Accelerometer 
Calibration  jJystem)  experiment  aboard  ar.  Agena  vehicle  in  a  near-polar  orbit. 

The  accelerometer  gave  density  values  with  an  estimated  accuracy  better  than  10Z 
and  a  time  resolution  of  J  s.  although  at  much  longer  and  irregular  intervals. 

■  ;  > 

The  dayside  response  (at  10.30  LT)  is  clearest  as  the  geomagnetic  index 
ap  starts  to  decline,  and  shows  peaks  on  the  magnetic  equator,  near  ±45°  geo¬ 
magnetic  latitude,  and  near  the  polar  caps.  There  are  pronounced  troughs  near 
±60°  geomagnetic  latitude;  in  the  northern  trough  (at  150km  height)  the  density 
falls  to  70Z  of  the  undisturbed  value.  In  succeeding  passes  the  structure 
collapses  irregularly.  The  equatorial  peak  nearly  disappears  within  3  h.  At 
the  peaks  the  density  in  near  its  maximum  as  ap  starts  to  decline,  indicating 
a  very  rapid  response;  but  the  density  in  the  troughs  increases  for  some  hours 
as  the  structure  collanses,  before  the  atmosphere  returns  to  normal. 

We  suggest  that  the  structure  shows  where  the  energy  is  injected,  at  least 
in  the  morning  sector,  in  an  intense  storm,  namely:  (i)  near  the  magnetic  equator, 
(ii)  near  the  particle  precipitation  Sones,  which  are  displaced  equatorwards  in 
the  storm,  (iii)  over  the  polar  caps.  Most  of  the  energy  input  is  probably  due 
to  particle  precipitation.  The  conclusion  of  DeVries,  that  most  of  the  heating 
occurs  in  two  bands  around  the  nightside  at  about  ±70°  geomagnetic  latitude,  and 
that  the  energy  is  transported  over  the  rest  of  the  earth  by  atmospheric  circula¬ 
tion  and  gravity  waves,  is  not  confirmed. 
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INTRODUCTION 


On  12  Kay  1967,  the  US  Air  Force  launched  into  a  cear-polar  orbit  an 
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attitude-controlled  Agena  vehicle  carrying  an  extremely  sensitive  accelerometer  * 
generally  known  as  LOGICS  (Lov-0  Accelerooeter  Calibration  System) .  Fortuitously, 
the  most  intense  geomagnetic  disturbance  since  I960  occurred  during  the  flight. 

Two  class  2B  solar  flares  and  one  class  3B  solar  flare  on  23  Kay  were  followed 
by  increased  soler  activity,  and  by  greatly  increased  geomagnetic  activity  on 
25-26  May.  As  Fig. I  shows,  the  I0.7ca  solar  radiation  flux  rose  steadily  from 
178.3  (in  units  of  10  ^  W  m  ^  s  * )  on  22  May  1967  to  213.2  on  26  May.  The 
3-hour ly  geomagnetic  index  a^  started  tx,  increase  after  revolution  41,  reaching 
a  maximum  of  236  on  revolution  45,  falling  to  154  on  revolution  47,  and  then 
increasing  again  to  400  on  revolutions  51,  52  and  (part  of)  53  before  slowly 
decaying  to  less  than  100  by  revolution  62.* 

The  orbit  was  near-polar  with  the  north-going  passes  on  the  nightside  at 
*v  22.30  LT,  and  the  south-going  passes  on  the  dayside  at  'v  10.30  LT.  The  argu¬ 
ment  of  perigee  was  initially  u  ~  140°,  falling  to  120°,  so  that  perigee  was 
attained  at  a,  40°  N  rising  to  60°  N  when  the  satellite  was  moving  southwards  on 
the  dayside  pass.  After  an  orbital  manoeuvre  during  revolution  18  on  23  May 
in  which  apogee  was  raised  by  'v  70  km,  the  orbit  had  an  inclination  of  91°. 5, 
and  perigee  and  apogee  heights  of  147  and  403  km  respectively,  which  decreased 
to  141  and  296  km  over  the  next  three  days.  The  lowest  point  on  the  orbit  rela¬ 
tive  to  the  real  Earth  (as  opposed  to  perigee  relative  to  a  spherical  Earth)  was 
144.7  km  at  *  32°  N  on  rev.  18  falling  to  137.7  km  at  'u  41°  N  on  rev.  65. 

The  LOGAC;-.  experiment  provided  values  of  atmospheric  density  with  a  time 

resolution  of  1  s,  although  at  tcuch  longer  intervals,  and  an  estimated  accuracy 

.  .  1  2 
of  vichfn  10Z.  These  values  have  been  used  by  DeVries  ’  to  study  the  neutral 

atmospheric  density  at  quiet  times,  and  sore  especially  its  response  to 

geomagnetic  activity. 

Increases  in  the  density  of  the  neutral  atmosphere  associated  with 

3 

geomagnetic  activity  were  first  noted  by  Jacchia  from  the  analysis  of  the  drag 
acceleration  of  Sputnik  3  rocket.  From  a  later  study  of  Injun  3,  Jacchia  and 


*  In  the  LOGACS  data  the  index  ap  is  given  as  varying  continuously  simply  by 
attributing  each  3h  value  to  the  centre  of  the  interval  and  making  a  linear 
interpolation,  as  shown  by  the  broken  curve  in  Fig.l.  The  use  of  the  resulting 
pseudo-a  index  has  been  continued  in  this  Report. 
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Slowey  reported  that  the  heating*  was  greater  by  a  factor  of  4  or  5  near  the 
auroral  senes  than  at  low  latitudes.  In  later,  more  extensive,  studies  using 
4  satellites  and  covering  80  events,  Jacchia  et  al reached  sore  guarded 
conclusions.  For  small  disturbances  (Kp  <  5)  the  heating  rate  was  15-25Z 
greater  for  |t|  >  55°  than  at  lower  latitudes,  while  for  larger  disturbances 
the  heating  rate  w«»  ocrr.e times  greatly  enhanced  in  the  auroral  zones.  The 
heating  appeared  to  be  a  little  larger  on  the  nightside  for  Kp  <  5  ,  with 
the  opposite  trend  for  larger  storms,  but  the  scatter  was  so  large  that  the 
difference  was  not  significant.  They  also  found  that  the  atmosphere  responded 
to  geomagnetic  activity  rather  morA  quickly  at  high  latitudes  with  a  delay 
t  -  5.8  ±  0.5  h  for  ({'{  >  55°  (av.  65°)  compared  with  a  delay  t  *  7.2  ±  0.3 h 
for  |e|  <  55°  (av.  25°). 

On  the  other  hand,  by  analysing  200  disturbances  in  the  altitude  range  250 
to  800  km  using  drag  data  from  6  satellites,  Roemer^  concluded  that  the  time 
delay  was  independent  of  latitude,  height  or  local  time;  on  average  the  ■mp™— «— 
density  occurred  5.5  ~  0.3  h  after  the  maximum  geomagnetic  disturbance.  He  also 
found,  however,  that  the  temperature  increase  was  greater  at  higher  latitudes, 
and  also  somewhat  greater  O'-  2QZ)  by  night  than  by  day. 

It  m^st  be  pointed  out  that  results  based  on  orbital  analysis  necessarily 
involve  averaging  over  a  few  orbits,  as  well  as  over  a  range  of  height  and 
latitude,  and  therefore  could  well  contain  an  intrinsic  time  delay.  In  the 
LOGACS  experiment,  it  was  possible  to  compare  the  in  situ  acceleration  data  with 
orbital  analysis  based  on  radar  observations  with  the  conclusion  that  the  atmos¬ 
pheric  response  deduced  from  orbital  analysis  has  an  inherent  delay  of  about  3  h. 
Although  the  authors'’’*’  might  not  agree,  it  could  well  be  that  the  previously 
quoted  values  of  time  delay  should  be  reduced  by  about  3  h. 

Against  this  somewhat  confused  background  of  results  on  a  presumably 
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complex  phenomenon,  DeVries  ’  came  to  the  following  conclusions  from  his 
analysis  of  the  LOGACS  data: 

(i)  During  quiet  times  the  density  is  lower  ('u  10Z)  than  expected  in  the  night¬ 
time  (north)  auroral  regions  and  higher  (y  30Z)  in  the  daytime  (north) 
auroral  regions. 


*  The  increased  density  appears  to,  or  is  presumed  to,  depend  on  height  in  much 
the  same  way  as  would  the  unperturbed  density  for  a  higher  level  of  solar  act¬ 
ivity  (as  measured  by  the  Fio.7  index,  which  is  assumed  to  be  an  indicator  of 
the  intensity  of  the  solar  EUV  flux);  the  change  in  density  can  then  be 
expressed  in  terms  of  an  equivalent  change  in  exospheric  temperature,  T*  , 
using  an  atmospheric  model. 


(ii)  The  density  in  the  night-tine  (north)  auroral  regions  increased  almost 

'  simultaneously  with  the  increase  in  geonagnetic  activity  -  certainly  with 
a  delay  less  than  90  m  (one  orbital  period).  At  the  sane  tine  the  density 
suffered  a  temporary  decrease  by  up  to  30 Z  in  the  daytime  (north)  auroral 
regions . 

(iii)  The  delay  time  vss  latitude-dependent,  increasing  from  nearly  simultaneous 
in  the  night-time  auroral  regions  to  from  4i  to  6  h  in  the  dayside 
equatorial  regions. 

DeVries  concluded  that  most  of  the  energy  provided  to  the  atmosphere 
during  enhanced  geomagnetic  activity  was  deposited  in  the  nightside  auroral 
regions,  in  fact  by  Joule  heating  at  or  below  150  km,  and  transported  from  there 
by  convective  circulation  and/or  atmospheric  waves.  Moreover  he  interpreted 
vavelike  patterns  of  density  variation,  with  wavelengths  of  the  order  of  450  km 
and  half-amplitudes  of  up  to  35Z,  as  evidence  in  favour  of  energy  transfer  by 
horizontally-propagating  waves.  The  aim  of  the  present  study  was  to  make  an 
independent  examination  of  the  LOGACS  data  and  in  particular  the  evidence  for 
atmospheric  waves. 

2  PREFERRED  NUMBERS 

In  order  to  organize  and  assess  the  data  (which  had  been  stored  on  magnetic 
tape)  we  attempted  to  remove  the  height-dependence  of  the  density  without 
recourse  to  a  detailed  atmospheric  model,  by  assuming  that  the  density  scale 
height  increased  linearly  with  height,  taking  the  density  and  the  scale  height 

-q  _3 

to  have  fixed  values  at  90  km,  namely  p^  *  3.46  x  iO  g  cm  and 

*  5.535  km.  The  validity  of  this  assumption  of  constant  dB/dz  can  be  judged 
from  Fig. 3  for  geomagnetically  quiet  times  above  140  km.  A  measured  density  p 
at  some  height  z  can  then  be  converted  to  a  density  p*  at  some  standard 
height  z*  ,  which  was  taken  in  fact  to  be  140  km  -  near  the  actual  height  of 
perigee.  The  values  of  the  reduced  density  p*  were  then  plotted  (by  automatic 
graph-plotting)  for  each  dayside  and  nightside  pass  separately.  There  was  in 
fact  a  dip  in  p*  around  the  actual  perigee  (on  dayside  passes)  indicating 
either  that  dB/dz  changes  significantly  belov  140  km,  or  that  the  fixed  values 
at  90  km  are  wide  of  the  mark.  The  giuphc,  however,  admirably  served  their 
purpose  of  assessing  the  data. 

The  graphs  contained  a  number  of  seemingly  artificial  features  which  were 
found  to  be  related  to  particular  values  of  density  in  the  original  data.  The 


qos t  striking  instance  occurred  on  the  nightside  pass  between  revs.  41  arid  42 
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at  heights  approaching  350  km.  Although  DeVries  *  would  assume  that  a'il 
measurements  above  300  ks  were  suspect,  it  will  serve  admirably  as  an  example. 

Thus  the  density  values  contained  the  following  sequence  of  numbers 
(all  *  10  ^  g  cm  ^): 

0.412,  (0.425),  0.412,  0.309,  0.258,  0.206,  (0.200),  0.206,  0.258,  0.309,  0.412, 
0.515,  (0.590),  0.515,  0.412,  0.309,  (0.307),  0.309,  0.412,  0.515,  (0.696), 

0.773  ....  . 

Clearly  the  same  set  of  numbers  (ignoring  those  in  brackets)  is  being  repeated. 

A  count  of  all  (roughly  6000)  density  values,  taking  the  3-digit  part  O.abc 

and  ignoring  the  exponent,  showed  that  the  numbers  0.103,  0.129,  0.155,  0.206, 

0.258,  0.309,  0.412,  0.515  and  0.773  (and  possibly  also  0.121,  0.122  and  0.123) 

appeared  much  too  frequently.  Since  there  might  be  a  'quantisation*  effect 

associated  with  measuring  very  low  densities  at  the  greatest  heights,  the  counting 

process  was  repeated  for  each  exponent  separately  with  the  values  given  in  the 
— N  -3 

form  O.abc  x  10  ‘  (g  cm  ).  Some  selected  parts  of  the  resulting  frequency 
count  are  shown  in  Table  1 . 

Table  I  -  Frequency  of  density  valuesf 


abc 

N  -  11 

N  *  '.2 

N  «  13 

N  «  14 

H  -  15 

102 

5 

6 

6 

0 

0 

103 

26 

2! _ 

46 

2 

0 

104 

10 

3 

3 

0 

s 

V 

120 

7 

4 

3 

0 

0 

121 

4 

34 

4 

0 

0 

122 

4 

50 

5 

0 

0 

123 

3 

44 

3 

1 

0 

124 

10 

9 

4 

0 

0 

128 

9 

3 

3 

0 

0 

129 

29 

29 

45 

3 

0 

130 

8 

4 

1 

0 

0 

"al 

f  Densities  are  given  in  the  fora  O.abc  x  10  ‘  g  cm 

The  conclusion  is  that  the  preferred  numbers  103,  129,  155,  206,  258,  309,  412 
515  and  773  are  associated  with  all  exponents;  those  underlined  are  multiples  of 
103,  and  the  reminder  are  related  to  them  as  nearly  as  possible  as  the  ratio  of 
small  integers.  For  exai?ole 

258  =  2  x  129  =  -|  *  103  =  |  *  515  =  j  *  773  ,  etc. 
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rhis  ccmirt  also  confirms  the  possible  triplet  already  noted  (namely  121 ,  122, 

123)  and  also  reveals  two  further  triplets  151,  152,  153  and  181,  182,  183;  all 

three  triplets  appear  such  too  frequently  but  only  in  association  with  one 
-1 1 

exponent,  naaely  10 

These  preferred  numbers,  18  altogether,  account  for  over  1000  of  the  roughly 
6000  density  values.  They  occur  aost  often  in  the  earlier  undisturbed  portion 
of  the  data  before  the  geoaagnetic  stora,  where  they  often  appear  as  quite 
'reasonable*  and  'realistic*  values;  the  effect  is  that  soae  of  the  points  show 
a  too  exac'-iy  regular  pattern  of  variation.  Since  the  preferred  numbers  also 
soaetiaes  produce  features  which  are  clearly  artificial,  we  have  chosen  to 
reject  the  'preferred  numbers  *  everywhere. 
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In  particular,  DeVries  '  has  put  forward  wavelike  patterns  of  density 
variation  as  evidence  in  favour  of  the  transfer  of  energy  fro*  the  auroral 
zones  by  horizontally-propagating  waves.  The  actual  patterns  occur  just  north 
of  the  ascending  nodes  (on  the  aightside)  on  revs.  53  and  54  at  heights 
between  1*50  and  300  km,  and  are  shown  in  Fig. 7  of  Ref. I  (revs.  53  and  54)  and  in 
Tif  .7  c  £  Fef.2  (rev.  54  only).  We  find  that  these  patterns  consist  largely  of 
'preferred  augers',  and,  in  our  opinion,  cannot  be  accepted  as  evidence. 

3  A  STAKIV-KD  OUIfcT-TDg  DENSITY 

vvjm  the  authentic  plots  of  the  'reduced  density*  of  section  2,  it  appeared 
that  the  most  complete  and  the  aost  structured  parts  of  the  data  were  to  be  found 
on  the  dayside  passes.  Accordingly  we  sought  a  more  satisfactory  way  of  quanti¬ 
fying  the  response  than  finding  a  'reduced  density*.  Bearing  in  mind  that  we  are 

m  S  ^ 

concerned  only  with  a  very  particular  section  through  '<he  atmosphere  defined  by 

local  time,  season  (possibly)  and  a  correlated  variation  of  height  with  latitude, 

it  ..'eeas  simplest  and  best  to  find  a  comparison  density  for  geomagneticallv  quiet 
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time:'  from  the  data  of  the  LOGACS  experiment  itself.  DeVries  *  in  fact  used  the 
densities  on  rev.  4!,  but,  since  the  characteristics  of  the  orbit  (local  time, 
perigee  pi  >it£on,  eccentricity)  change  little  over  the  four  days,  it  is  possible 
to  combine  all  the  data  from  the  first  41  orbits.  Fig. 2  shows  all  1718  values 
of  log  p  plotted  against  height  from  the  dayside  passes  (defined  by  local 
time  lying  between  S.9  and  11.5  b)  up  to  revolution  41,  with  only  the  'preferred 
numbers*  omitted.  Clearly  the  density  is  a  'good*  function  of  height  only,  with 
little  scatter  below  about  250  ks,  although  the  passes  include  some  with  slight 
to  moderate  geomagnetic  activity  (a  maximum  of  a  >31  on  rev.  33). 
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The  lowest  heights  from  140  to  170  km  refer  to  the  region  around  perigee 

fro®  the  equator  to  the  (north)  daytime  auroral  zones,  while  the  greatest  heights 

are  attained  approaching  the  South  Pole.  An  inspection  of  the  plots  of  ’reduced 

density*  shows  that  the  scatter  above  250  km  arises  because  the  density  pole- 

wards  of  about  40°  S  is  extremely  variable  from  pass  to  pass.  For  example  the 

’reduced,  density’  approaching  the  South  Pole  shoots  up  on  rev.  34  and  falls  to 

nearly  zero  on  rev.  35.  This  variability  does  not  see®  to  be  correlated  with  any 

variations  in  a  .  Indeed  it  is  possible  that  these  features  should  be  treated 
P 

as  fictitious.  Certainly  the  density  patterns  often  seem  too  regular  for  ecaforfc, 
I  2 

and  DeVries  *  would  reject  all  measurements  above  300  km  as  suspect.  On  the 
other  hand  this  variability  is  not  inconsistent  with  the  results  of  Blaaont  and 
Luton^,  who  find  that  the  temperature  at  F-region  heights  shows  unexpected  maxima 
in  the  polar  regions  even  at  geomagnetical ly  quiet  times. 

Ail  the  daytime  values  of  log  p  up  to  rev.  41  (with  only  the  ’preferred 

nutters'  omitted)  as  shown  in  Fig. 2,  are  well-fitted  by  a  cubic  polynomial  in 

the  height;  the  sum  of  squares  of  the  residuals  is  only  marginally  reduced  in 

going  beyond  a  cubic.  Scsse  of  the  individual  quiet  daytime  passes  have  been 

compared  with  the  density  generated  from  the  cubic  polynomial.  In  fact  rev.  41, 
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which  was  used  by  DeVries  is  the  most  nearly  complete  pass,  and  also  corres¬ 
ponds  very  closely  with  the  polynomial  fit.  Fig. 3  shows  the  values  of  the 
density  pseudo-scale  height  H  ,  with  error  bars  corresponding  to  ±lc  ,  derived 
from  the  cubic  fit.  The  value  of  E  increases,  at  first  nearly  linearly,  from 
20  km  at  140km  height  to  a  maximum  of  about  38  km  at  270km  height  and  thereafter 
falls  again.  Since  the  density  values  refer  to  a  particular  section  through  the 
atmosphere,  H  is  not  the  scale  height  at  a  given  place.  Thus  it  is  not  being 
suggested  that  the  density  scale  height  derived  from  a  vertical  section  through 
the  South  polar  regions  is  less  than  30  ks  at  350to  height;  instead  it  will  lie 
somewhere  between  the  value  30  to  and  an  extrapolation  of  the  early  linear  part 
of  the  curve  in  Fig. 3.  The  linear  part  of  the  curve  in  Fig. 3  is  probably  a  good 
measure  of  the  (vertical)  density  scale  height  up  to  mid-latitudes.  Nevertheless, 
if  the  observations  above  250  to  can  be  accepted.  Fig. 3  indicates  that  the  atmos¬ 
phere  in  the  South  polar  regions  in  May  is  on  balance  colder  (although  on  sum 
passes  it  is  hotter) . 

4  RESPONSE  OF  THE  DAYTIME  ATMOSPHERE 

The  greatest  and  the  most  clearly  developed  response  of  the  daytime  atmos¬ 
phere  occurs  just  after  the  main  peak  of  the  geomagnetic  storm  when  reaches 


400  for  3  h  (rev.  51  Co  rev.  53).  Fortunately  the  daytime  pass  on  rev.  53  gives 
a  nearly  complete  coverage  in  latitude,  and  rev.  55  is  even  more  complete.  The 
response  is  shown  in  Figs. 4  and  5  in  the  form  of  the  ratio  p/Pq  »  where  p^  is 
the  standard  quiet-time  density  derived  in  section  3,  plotted  against  geographic 
latitude.  Fig. 4  refers  to  the  daytime  pass  on  rev.  53,  during  which  a^  falls 
from  400  to  383;  there  is  a  gap  in  the  latitude  coverage  from  50°  S  to  80°  S. 

Fig. 5  refers  to  rev.  55  during  which  a^  falls  from  306  to  276,  and  is 
essentially  complete. 

The  two  figures  are  very  similar;  excepting  the  gap  in  Fig. 4,  they  show  a 
wave-like  pattern  with  a  central  maximum  at  10°  N,  further  maxima  at  50°  K  and 
40°  S,  and  yet  further  maxima  approaching  the  polar  caps.  The  pattern  is  not 
symmetrical  in  geographic  latitude,  but  is  much  more  cearly  symmetrical  (in  form 
if  not  in  magnitude)  in  magnetic  coordinates;  the  magnetic  equator  and  the  feet 
of  various  L-shells  are  marked  in  Figs. 4  and  5.  In  magnetic  coordinates  the 
central  peak  occurs  nearly  at  the  equator,  and  the  next  peaks  are  near  L  =  2  in 
the  northern  hemisphere,  but  perhaps  nearer  L  -  3  in  the  southern  hemisphere. 
This  slight  asysaetry  between  north  and  south  in  the  location  of  the  peaks  and 
tro-  hs  could  be  due  to  the  greater  height  in  the  southern  hemisphere;  alter¬ 
natively  the  geometry  of  the  situation  could  be  imposing  a  compromise  between 
* 

geomagnetic  and  geographic  latitudes.  On  the  other  hand  there  is  little  doubt 
that  the  greater  magnitude  of  p/Pq  in  the  southern  hemisphere  is  simply  due  to 
the  greater  height. 

Although  the  gross  features  of  the  response  have  a  wave-like  pattern,  the 
gross  structure  is  not  a  propagating  wave.  The  dayside  pass  on  rev. 54  is  some¬ 
what  incomplete  in  coverage,  but  what  there  is  lies  neatly  intermediate  between 
the  patterns  of  Figs. 4  and  5.  In  the  succeeding  dayside  passes,  the  structure 
collapses  in  irregular  fashion.  It  is  particularly  noticeable  that  the  equatorial 
peak  has  already  nearly  disappeared  by  rev.  55. 

There  are  certainly  small-scale  wave-like  features  such  as  the  double  peak 
at  50°  N  in  Ifig.5,  but  we  would  hesitate  to  interpret  them  as  evidence  for 
propagating  waves.  Unfortunately  there  is  not  consistently  enough  data  with 
sufficiently  close  coverage  to  be  able  to  identify  small-scale  features  on 
successive  passes.  In  any  case  it  must  he  remembered  that  each  successive  pass 
samples  the  atmosphere  about  23°  further  vest. 


We  suggest  that  the  structure  of  Figs. 4  end  5  shot1-  si o-”  v  “he  Jfr.'] 

is  injected  into  the  atmosphere,  at  least  in  u.  **  •  .i.u~.ig  sector,  in  5-  '-ns.- 

geomagnetic  event,  namely: 

(i)  Near  the  magnetic  equator. 

(ii)  In  the  particle  precipitation  zones  which  are  displaced  towards  tic 
equator  in  a  geomagnetic  store.  The  sore  intense  the  tom,  the  greater 
is  the  displaceeent.  As  already  noted,  this  peak  app  k*S  near  L  *  2  in 
the  northern  hemisphere  but  nearer  L  *  3  in  the  southern  J-rsn sphere. 

(iii)  Over  the  polar  caps. 

One  of  the  features  of  the  atmospheric  response  is  the  appearanc  of 
pronounced  troughs  at  about  ±60°  magnetic  latitude,  i.e.  approaching  nor. ^-1 
auroral  latitudes.  The  northern  trough  appears  in  both  Figs. 4  and  5,  and  rev 
a  minimum  density  of  only  70Z  of  the  quiet-time  value  at  a  height  of  !50  km.  b 
corresponding  trough  in  the  southern  hemisphere  is  missing  on  rev.  53  doe  to 
the  gap  in  coverage,  but  shows  clearly  in  Fig. 5.  Here  the  sisiam  value  of  p/p, 
is  about  1.3  at  a  height  of  270  km.  Two  possible  reasons  for  these  troughs  axe 
as  follows: 

(i)  The  air  has  been  accelerated  and  displaced  by  electrodynamic  forces. 

(ii)  The  air  is  colder,  due  to  the  removal  of  part  of  its  usual  heat  source  by 
the  displacement  of  the  particle  trapping  boundaries  to  much  lower 
latitudes. 

In  either  case  it  is  assumed  that  the  trough  at  the  greater  height  in  the  soothe 

hemisphere  has  been  ‘f tiled -in '  to  some  extent  by  spill-over  from  the  adjacent 

peaks.  On  the  other  hand  the  densities  in  the  troughs  appear  to  be  in  agreement 

8 

with  the  theoretical  predictions  of  Mayr  and  Vo Hand  .  Assuming  that  heating 
takes  place  near  the  base  of  the  thermosphere  (say  between  100  and  150  km),  they 
find  that  the  thermospheric  wind  circulation  excited  during  magnetic  storms 
removes  atomic  oxygen  at  high  latitudes.  The  wind-induced  variations  in  atomic 
oxygen  exceed  temperature  effects  below  250  km  giving  a  reduction  in  density. 

The  calculated  depletion  is  most  pronounced  near  180  km  where  the  density  can 
decrease  by  a  factor  of  2.  Above  250  ka  the  rise  in  temperature  predominates 
and  the  density  increases. 
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5  POSSIBLE  HEATING  MECHANISMS 
2 

DeVries  discussed  the  possible  means  by  which  energy  is  deposited  in  the 
thermosphere  during  geomagnetic  disturbances  in  terms  of  a  i -dimensional  treat- 

Cl 

ment  by  yolland'  in  which  horizontal  transport  of  mass,  momentum  and  energy  are 
neglected.  Volland  considered  three  possible  locations  for  the  original  heating: 

(i)  Heating  at  the  base  of  the  thermosphere  by  Joule  dissipation  (or  wave 

IQ 

energy).  Cole  \  who  has  particularly  emphasised  the  possible  effects  of 
Jsule  heating  in  the  auroral  zones,  estimated  that  the  maximum  heating 
takes  place  near  150  km. 

(ii)  Energy  input  throughout  the  entire  thermosphere.  This  would  correspond, 
in  particular,  to  heating  by  particle  precipitation  as  suggested  by 
Mcllwain^,  which  would  give  a  more  or  less  uniform  rise  in  temperature 
(assuming  the  cross-section  is  independent  of  velocity)  down  to  the  level 
at  which  the  particles  are  stopped. 

(iii)  Heating  at  the  top  of  the  thermosphere,  which  would  correspond  to,  for 

instance,  heating  by  hydromagnetic  waves  generated  by  the  interaction  of 

1 2 

the  solar  wind  with  the  magnetosphere  as  suggested  by  Dessler  . 

9 

Volland  concluded  that  the  geomagnetic  activity  effect  was  due  to  particle 
precipitation  into  the  nightside  auroral  zones,  i.e.  heat  source  (ii),  with  a 
subsequent  distribution  in  latitude  and  longitude  by  conduction  and  atmospheric 
circulation.  This  would  be  in  general  accord  with  the  conclusions  from  satellite 
drag  observations  of  an  enhanced  response^  and  a  shorter  time  delay"*  at  high 
latitudes. 

Since  he  found  an  immediate  increase  in  density  over  the  wide  latitude 
band  polewards  of  about  50°  N  on  the  nightside,  whereas  there  was  a  decrease  at 
70°  N  on  the  dayside,  DeVries  reached  the  same  conclusion  on  the  location  of 
the  original  heating.  He  attributed  the  heating,  however,  to  Volland's  heat 
source  (i)  rather  than  (ii);  to  be  specific,  DeVries  concluded  that  the  mechan¬ 
ism  was  Joule  dissipation  (at  about  150  km)  around  the  night-time  auroral  zones, 
with  the  energy  being  transported  to  the  rest  of  the  earth  by  atmospheric  circula- 
tion,  conduction,  and  gravity  waves,  with  the  emphasis  perhaps  on  this  last 
mechanism.  He  noted,  however,  that  an  immediate  response  over  a  wide  latitude 
band  on  the  nightside.  (at  about  200ko  height)  is  perhaps  at  variance  with  the 
theoretical  predictions  of  Chiu  and  Ching  ,  In  studying  the  effectiveness  of 
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gravity  waves  in  distributing  an  energy  input,  e.g.  from  Joule  heating,  they 
concluded  that  there  was  a  density  time  lag  of  6  ±  2  h  even  directly  above  the 
heat  source. 

The  present  results  on  the  dayside  response  in  the  morning  sector  do  not 

2 

support  the  conclusions  of  DeVries  .  As  already  noted,  the  response  has  a  well- 

developed  structure  by  rev.  53  (when  a  is  just  beginning  to  decrease  from  the 

P 

maximum  value  of  400  attained  3  h  earlier)  with  a  peak  near  the  magnetic  equator, 
further  peaks  at  about  ±45°  magnetic  latitude,  and  yet  further  peaks  approaching 
the  polar  caps.  The  response  at  the  peaks  appears  to  be  near  itr  maximum  on 
rev.  53,  so  that  the  time  delay  cannot  be  long.  On  the  other  hand  the  density 
in  the  troughs  continues  to  increase  for  some  hours  as  the  pattern  collapses 
and  the  troughs  are  filled  in,  before  the  atmosphere  eventually  begins  to  return 
towards  normality.  Since  the  peak  values  of  p/Pq  in  the  southern  hemis  Here 
exceed  the  corresponding  peaks  in  the  northern  hemisphere  (which  are  at  loter 
heights),  the  implication  would  appear  to  be  that  a  substantial  part  of  the 
energy  input  must  raise  the  temperature  fairly  uniformly  throughout  the  height 
range,  in  other  words  that  the  physical  mechanism  is  probably  particle 
precipitation, 

6  CONCLUSIONS 

The  original  aim  of  the  present  study  had  been  to  make  an  independent 
examination  of  the  LOGACS  data,  and  in  particular  of  the  possible  role  of  gravity 
waves  as  a  means  of  transporting  energy.  In  the  event,  however,  our  investiga¬ 
tions  have  taken  a  rather  different  course.  Our  initial  survey  of  all  the  data, 
taking  all  dayside  (  m  10.30  LT)  and  nightside  (  *  22.30  LT)  passes  separately, 
showed  that  the  most  complete  and  the  best  ordered  data  were  to  be  found  on  the 
dayside  passes.  The  dayside  passes  also  have  the  advantage  of  showing  the 
response  in  both  hemispheres;  this  is  not  possible  for  the  nightside  passes 
since  the  satellite  is  above  300  km,  where  the  data  is  sparse  and  unreliable,  in 
the  southern  hemisphere. 

The  most  clearly  developed  dayside  response  occurs  just  after  the  main 
peak  of  the  storm.  Figs. 4  and  5  show  the  density  ratio  p/Pq  »  where  Pq  is  a 
standard  quiet-time  density  derived  from  the  earlier  part  of  the  LOGACS  data, 
for  rev.  53  (just  as  a^  starts  to  decline  from  its  peak  value  of  400)  and 
rev.  55  respectively.  Rev.  54  is  not  shown  because  it  is  incomplete.  The 
response  is  nearly  synmetrical,  in  form  if  not  in  magnitude,  in  terms  of 
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magnetic  coordinates,  with  a  peak  near  the  magnetic  equator,  further  peaks  near 
±45°  magnetic  latitude,  and  yet  further  peaks  approaching  the  polar  caps.  There 
are  also  pronounced  troughs  near  ±60°  magnetic  latitude.  The  northern  trough 
(at  150km  height)  appears  in  both  figures;  on  rev.  53  the  density  is  only  702  of 
the  undisturbed  value,  but  has  recovered  to  852  by  rev.  55.  The  corresponding 
southern  trough  is  missing  on  rev.  53  due  to  a  gap  in  coverage,  but  appears 
clearly  on  rev.  55  where  the  density  (at  a  height  of  270  km)  has  Gnly  increased 
by  a  factor  of  1.3  from  its  undisturbed  value.  At  the  peaks,  the  density  is 
near  its  maximum  on  rev.  53  indicating  that  the  time  delay  in  the  response  is 
relatively  short.  In  particular  the  equatorial  peak  on  rev.  53  has  nearly 
disappeared  three  hears  later,  perhaps  because  Ion  drag  offers  negligible  resist¬ 
ance  to  an  expansion  in  latitude.  On  the  other  hand  the  density  in  the  troughs 
continues  to  increase  for  some  hours;  in  succeeding  passes  the  structure  collapses 
in  irregular  fashion  -  indeed  'disintegrate*  might  be  a  more  apt  verb. 

We  suggest  that  the  dayside  response  simply  shows  where  the  energy  has 
been  injected  into  the  atmosphere,  at  least  in  the  morning  sector,  during  an 
intense  storm,  namely j 

(i)  near  the  magnetic  equator, 

(ii)  near  the  particle  trapping  boundaries  which  are  displaced  equ&torwards  in 
a  geomagnetic  storm, 

(iii)  over  the  polar  caps. 

In  a  further  paper  by  one  of  us,  it  will  be  argued  that  the  equatorial  response 

is  related  to  the  decay  of  the  ring  current.  This  conclusion  is  in  agreement 

14 

with  the  results  of  May  and  Hiller  on  the  response  of  the  near-equatorial 

atmosphere  (at  *v  310  km  near  noon  local  time)  from  the  decay  of  the  spin-rate 

of  Ariel  2  in  March-April  1964.  They  found  that  the  fine  structure  cf  the 

density  variations  was  more  closely  correlated  with  the  index  D  ,  an  hourly 

St 

index  based  only  on  the  horizontal  H-components  at  low-to-middle  latitude 
stations  and  which  el''  ctively  measures  the  strength  of  the  ring  current,  chan 
with  a 

P 

Turning  to  ~ht  mid-latitude  peaks,  the  density  ratio  £>/Pq  is  greater  at 
the  greater  heights  in  the  southern  hemisphere  and  the  peak  appears  just  as 
quickly.  The  conclusion  would  appear  to  be  that  the  temperature  must  increase 
fairly  uniformly  throughout  the  height  range,  and  that  a  substantial  part  of  the 
energy  input  is  probably  supplied  by  particle  precipitation. 
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It  is  not  surprising  that  orbital  drag  analysis  has  led  to  confusing  and 
contradictory  conclusions  if  the  complexity  of  the  day  side  response  is  typical. 
No  doubt  the  latitude  dependence  of  the  response  varies  with  local  tine  since  it 
is  unlikely  that  the  pattern  of  currents  or  the  pattern  of  particle  precipitation 
will  be  axially  syue-stric.  Also  the  response  might  reasonably  be  expected  to 
depend  on,  for  instance, 

(i)  the  inclination  of  the  magnetic  dipole  axis  to  the  earth-sun  line,  giving 
a  dependence  on  UT/longitude, 

(ii)  the  sense  of  the  interplanetary  magnetic  field,  including  its  north-south 
component, 

(iii)  the  magnitude  of  the  storm. 


Recently  further  in  situ  results,  which  all  refer  to  daytime  passes  near 

16.00  LX,  have  become  available  from  two  experiments  on  the  near-polar  satellite 

IS 

0G0  6.  Taeusch  et  at .  "  have  studied  the  response  during  several  magnetic 
storms  (with  a  maximum  of  80)  at  the  end  of  September  1969,  using  measure¬ 
ments  from  the  neutral  mass  spectrometer.  They  found  that  the  increase  in 
density  Cat  or  above  400  km)  was  more  or  less  confined  to  high  latitudes  with  the 
response  increasing  with  latitude,  and  the  atmosphere  being  scarcely  affected 
near  the  equator.  The  response  was  symmetric  in  terms  of  magnetic  latitude,  and 
the  time  lag  was  less  than  lh  polewards  of  ±50°.  Blamont  and  Luton7  also  report 
observations  covering  the  same  period  from  the  Fabxy-Perot  interferometer.  Here 
the  neutral  temperature  at  roughly  F-region  heights  is  derived  from  the  Doppler 

a 

broadening  of  the  6300  A  oxygen  line  in  the  day-glow.  They  find  that  the  res¬ 
ponse  to  geomagnetic  activity  is  strongest  in  the  polar  regions,  with  a  maximum 
temperature  of  about  1/00  K,  and  the  temperature  increase  falls  off  smoothly  to 
nearly  zero  at  the  equator.  They  also  report  unexpected  temperature  maxima  in 
the  polar  regions  even  at  geomagnetical ly  quiet  times. 

In  the  present  study,  where  the  LGGACS  data  before  the  storm  was  used  tc 
derive  a  comparison  undisturbed  density,  the  quiet-tin*  density  polewards  of 
40°  S  (and  above  25C  ka)  is  quite  variable  from  pass  to  pass  and  the  variation 
does  not  appear  to  be  correlated  with  the  (slight)  geomagnetic  activity  taking 
place,  at  least  as  measured  by  the  index.  This  variability  accounts  for 

the  scatter  above  250  ks  evident  in  Fig. 2,  and  is  not  inconsistent  with  the 
findings  of  Blamont  and  Luton7 .  If  the  observations  above  250  km  can  be 
accepted,  the  atmosphere  in  the  south  polar  regions  in  Hay  is  on  balance  colder 
chan  the  rest  of  the  days ids  thermosphere. 


Finally,  the  present  study  has  shown  that  the  LOGACS  data  is  imperfect  in* 
that  about  ?QG0  of  the  6QG0  density  values  are  accounted  for  by  18  'preferred 
numbers'  out  of  a  total  or  900  possible  3-digit  nuabers.  The  details,  given 
sore  fully  in  section  2,  are  quite  intriguing: 

(i)  9  of  the  numbers  are  related  as  nearly  as  possible  as  the  ratio  of  small 
integers,  and  appear  with  all  possible  exponents, 

(ii)  the  regaining  9  numbers  sake  up  3  sets  of  triplets  (121,122,123;  151,152, 
153;  181,182,183}  and  appear  only  with  one  exponent. 

Ho  explanation  has  been  discovered.  These  'preferred  nuabers*  occur  most  often 

in  the  earlier  undisturbed  part  of  the  data,  where  they  usually  appear  as 

'reasonable'  and  'realistic'  values,  but  shew  a  too  exactly  regular  pattern  of 

variation.  These  'preferred  numbers'  also  give  obviously  artificial  features, 

however ,  and  we  have  rejected  then  completely  in  our  analysis.  In  particular 

the  wave-like  patterns  of  density  variation  just  north  of  the  equator  on  the 

1  2 

nightsi.de  passes  on  revs.  53  and  54  which  were  put  forward  by  DeVries  *  et 
evidence  for  horizontally-propagating  waves,  consist  largely  of  these  ’preferred 
nuabers'  and,  in  our  view,  cannot  be  accepted.  It  is  possible  that  further 
imperfections  still  remain  in  the  data  -  in  places  there  is  still  a  tendency 
towards  a  too  exactly  regular  pattern  of  variation  -  but  this  has  not  been 
proved.  Since  the  LOGACS  experiment  has  provided  a  unique  opportunity  to  study 
the  atmospheric  response  to  an  intense  storm,  it  would  be  highly  desirable  to 
repeat  the  processing  of  the  original  raw  data  if  possible. 
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Geomagnetic  planetary  3-hour  index  qp,  and  mean  daily  solar 
radio  flux  at  IO-7cm,  during  Logacs  experiment 
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Fig.  2  Scatter  plot  of  quiet-time  Wayside  densities 
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Fig, 3  Density  pseudo-scale  height  H  from  quiet 


